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Abstract 

Conditions for the existence and stability of de Sitter space in modified gra- 
vity are derived by considering inhomogeneous perturbations in a gauge-invariant 
formalism. The stability condition coincides with the corresponding condition for 
stability with respect to homogeneous perturbations, while this is not the case in 
scalar-tensor gravity. The stability criterion is applied to various modified gravity 
models of the early and the present universe. 
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1 Introduction 



The 1998 discovery that the expansion of the universe is accelerated, obtained with 
the study of type la supernovae jT], has prompted many theoretical models to explain 
this phenomenon. Most of these models can be classified in three classes: dark energy 
models, modified gravity models, and brane-world models. In the first class it is assumed 
that a form of dark energy or quintessence of unknown nature has come to dominate the 
dynamics of the universe at recent times (redshifts z < 1). These models are usually 
explored in the context of Einstein's theory of general relativity [2|, or possibly in its 
scalar-tensor generalizations (extended quintessence) [3]. Dark energy must necessarily 
have exotic properties in order to generate acceleration. In the spatially flat Friedmann- 
Lemaitre-Robertson- Walker ("FLRW") line element describing our universe according 
to the recent cosmic microwave background experiments JI|, and given by 

ds 2 = -dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) (1.1) 

in comoving coordinates (t,x,y,z), the acceleration equation [5] 

- = > + 3P) (1.2) 
a b 

holds, where k = 8ttG and p and P are the total energy density and pressure of the cosmic 
fluid, respectively. An overdot denotes differentiation with respect to the comoving time 
t. If dark energy were the only form of energy of the universe, acceleration a > would 
require an exotic negative pressure P < — p/3; and if ordinary matter and dark matter 
contributing to the dynamics are taken into account, the pressure of dark energy must be 
even more negative to compensate. As a matter of fact, dark energy is even more exotic: 
the best fit to the supernovae data favours an extreme form of dark energy called phantom 
energy or superquintessence with P < —p - - or with an effective equation of state 
parameter w = P/p < —1. Furthermore, the effective equation of state should evolve 
with time [7j. If confirmed, this fact would definitely rule out the cosmological constant 
A as an explanation of the cosmic acceleration because w\ = — 1 is strictly constant (the 
cosmological constant model is anyway disfavoured because of the cosmological constant 
problem [HJ and of the cosmic coincidence problem |5| that accompany it). Most models 
of dynamical dark energy are based on a scalar field rolling in a potential V(<p), a way 
to implement cosmic acceleration that is well known from inflationary theory in the early 
universe jTOj- However, a canonical scalar field minimally coupled to the Ricci curvature 
R in Einstein gravity cannot explain an equation of state parameter w < — 1, which 
is equivalent to superacceleration H > 0, where H = a/ a is the Hubble parameter. 
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(This name distinguishes a regime in which the Hubble parameter increases from an 
"ordinary" acceleration regime in which a = a + H 2 ^j > and H < 0.) In fact, the 
energy density and pressure of such a scalar field are 



(.2 ±2 



p* = Y + v W> p + = y- v M> (L3) 

and the Einstein-Friedmann equation of general relativity 

H = -- Q {P + p) = -^ (1.4) 

yields H < for a universe dominated by such a scalar (the upper bound H =const. is 
attained by de Sitter space). In order to model an equation of state parameter w < — 1 
corresponding to H > 0, which is the situation favoured by the observational data, a 
phantom field with negative kinetic energy jTTj or a scalar field coupled non-minimally 
to gravity [T2] have been used. Both of these theories can be seen as special cases of 
scalar-tensor gravity, described by the action 



S = d^xyf^g 



- ^ s ab v a 0V fe - v(4>) 



1.5) 



where ip(4>) and u(<f>) are arbitrary coupling functions. The exotic properties of dark 
energy or of its extreme form, phantom energy, have led some authors to a different 
approach and to the second class of models mentioned above. Instead of postulating 
an exotic form of dark energy of mysterious nature, these authors consider the 

possibility that gravity deviates from Einstein gravity at large scales, and assume that 
the Einstein-Hilbert Lagrangian is modified by corrections that become important only 
at late times in the history of the universe, i.e., when the curvature becomes small. This 
class of theories, called modified gravity, is described by the gravitational action 



S 9 = j d* Xv ^g-f(R), (1.6) 

where f(R) is a non-linear function of R. The first model proposed had the form 
f(R) = R — fi A /R, in which the correction in R~ Y becomes important only at low 
curvatures R — > 0. The general form (|1.6J) of the action also includes quantum gravity 
corrections to Einstein's theory originally introduced to improve renormalizability fSJ 
|T5| and used in inflationary models of the early universe In addition to the desired 
phenomenological properties of modified gravity in cosmology, there is some motivation 
for these models from M-theory [TH] . 
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In both classes of models, depending on the arbitrary functions and parameters 
adopted, there are solutions describing universes that accelerate forever, other solutions 
in which the universe ends its existence in a finite time in the future in a Big Rip 
singularity |19 3 or encounters another type of "sudden future singularity" [2U]. The fate 
of the universe depends on whether attractor solutions that are forever accelerating, or 
Big Rip attractors exist in the phase space, and on the size of their respective attraction 
basins. In many models of both dark energy and modified gravity there are de Sitter 
attractors accelerating forever. In this paper we focus on modified gravity and in sec. 3 
we use scalar-tensor gravity for a comparison of properties. We determine whether de 
Sitter attractors exist in the phase space of modified gravity by deriving conditions for 
the existence and stability of these solutions. Throughout most of this paper we consider 
general non-linear actions of the form with d 2 f/dR 2 7^ and, in the final part of 
this paper, we apply our general results to specific scenarios and forms of f(R) proposed 
in the literature. 

While it is straightforward to study the stability of de Sitter space with respect 
to homogeneous perturbations, which only depend on time, it is physicallly more sig- 
nificant to assess stability with respect to more general inhomogeneous perturbations, 
which depend on both space and time. This goal is more ambitious because of the 
gauge-dependence problems associated with this type of cosmological perturbations [21] 
and one expects the stability condition with respect to inhomogeneous perturbations to 
be more restrictive than the corresponding condition for stability with respect to homo- 
geneous perturbations. It comes therefore as a surprise that these conditions coincide, 
as shown in sec. 3 and briefly reported in a previous communication [22] _ this result 
can not be guessed or justified a priori. A similar analysis shows instead that in scalar- 
tensor theories the stability condition with respect to homogeneous perturbations is 
indeed more restrictive than the corresponding one for homogeneous perturbations [22\. 

Certain modified gravity models are ruled out on the basis of instabilities that mani- 
fest on short timescales (221 123 1213 EHl EH] • The stability condition derived here has the 
advantage of being applicable to any non-linear Lagrangian of the form \f—g f{R) and 
is useful in the study of the phase space and dynamics of modified gravity scenarios. 
Another motivation for our study is that, in order to be viable, modified gravity models 
need to have the correct Newtonian and post-Newtonian limit, and currently there is 
disagreement on whether certain models pass or not this test |2Hl 123 I2H] • Because 
many models do not admit a Minkowski solution around which to expand the weak-field 
metric, an expansion around the de Sitter background is used instead [211123 EDI- This 
is meaningful when a de Sitter solution exists and is stable. 

The plan of this paper is the following: in sec. 2 we derive a stabiliy condition for 
de Sitter space with respect to inhomogeneous perturbations by using a covariant and 
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gauge-invariant formalism suitable for generalized gravity (including scalar-tensor and 
modified gravity, and possibly mixed models). In sec. 3 we derive the much simpler 
stability condition with respect to homogeneous perturbations in modified gravity, and 
we interpret our results. In sec. 4 the stability condition derived for modified gravity 
is applied to various scenarios widely discussed in the literature, while sec. 5 contains a 
discussion and the conclusions. 

2 Stability of de Sitter space with respect to inho- 
mogeneous perturbations 

We can consider at once modified gravity and scalar-tensor theories, which we will use 
for a comparison with modified gravity, by studying the gravitational action 



S = d A x \f—g 



(2.1) 



This action contains also possible combinations of modified and scalar-tensor gravity if 
both Of /dip and d 2 f/dR 2 are non-vanishing: such mixed scenarios have received little 
attention in the literature so far [32]. Scalar-tensor gravity is the special case in which 
/ is linear in R, i.e., f((f>,R) = tj)((j))R, while modified gravity corresponds to setting 
= 1 and d 2 f/dR 2 ^ 0. In the spatially flat FLRW metric (jl.lj) the field equations 
assume the form 

H 2 = —(^<p 2 + ^-L + V-3HF) , (2.2) 
3F V 2 2 2 / ' K J 

H = -^(u^ + F-HF) , (2.3) 
1 fduj - 2 df , n dV\ 



* +3 ^ + ^U^ 9T 2 ¥r ' (2A) 

where F = df /OR. It is natural to use (H, 0) as dynamical variables and the equilibrium 
points of the dynamical system ()2.2|) - (j2.4j) are de Sitter spaces with constant scalar field 
(Ho, 0o): they exist subject to the conditions 

#0^0 = 2 (/o - V ) , (2.5) 
fo = 2^' , (2-6) 
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where Rn 



12#o > /o 



y, = dV 
V — dcj> 



and a 



f(<h,Bo), F = F (0o, -Ro), Vq = V(0 O ) 

prime denotes differentiation with respect to <fi. In modified gravity, there is only the 
condition ()2.5|) for the existence of de Sitter solutions because there is only one arbitrary 
function f(R). 

In order to describe inhomogeneous perturbations of the de Sitter fixed points (H , O ) 
we use the covariant and gauge-invariant formalism of Bardeen-Ellis-Bruni-Hwang- 
Vishniac [2U ESj in the version given by Hwang and Hwang and Noh [31] for generalized 
gravity. The metric perturbations A, B, Hi, and Ht are defined by the relations 



#00 



-a 2 (1 



2AY) 



(2.7) 



9oi 



-a 2 B Yj 



{21 



9ij 



a 2 [h ij {l + 2H L )+2H T Y i 



(2.9) 



where the scalar harmonics Y are the eigenfunctions of the eigenvalue problem Vj V* Y = 
—k 2 Y, and where hij is the three-dimensional metric of the FLRW background, Vj is 
the covariant derivative associated with h^, while k is the eigenvalue. The vector and 
tensor harmonics Yi and satisfy the equations 



Y, 



1 



k 

1 _ _^ 



1 



Yij = ■^V i V J Y + -Yh 
We need the Bardeen gauge-invariant potentials 



-(aH T 



and the Ellis-Bruni variable jHSj 

A0 = 



(2.10) 
(2.11) 



(2.12) 
(2.13) 
(2.14) 



(2.15) 
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while equations similar to eq. (|2.15|) define the gauge-invariant variables Af, AF, and 
AR. The first order equations satisfied by the gauge-invariant perturbations are 



UJ 



k 2 ft d (I du\ d f 1 df 1 dV\ 

a 2 2 d(p \u d<p ) d<p\2oj dcj) uj d(j) J 



- 3$ 



*A(df dv\ i d 2 f 

uj \d<f) d<f>) 2u d<pdR 



A0 
(2.16) 



k 2 R 



1 f l2 du df A dV\ 



AF + 3HAF +(— — — ) AF + — AR + - cu0A0 + - ( ft — + 2^-4 A0 

(2.17) 



F[$ A -3$ H ) + -(FR-2f + AV) $ A 



H T + \3H + - \H T + — H T = 



a" 



(2.18) 



F \ 1 f AF AF u ■ 

®H+\H + —j<S> A = - ^— - H— + - A0 



(2.19) 



k\\ 1 / uj - 2 3F 2 \ 1 ]3FAF / • A; 2 3H F \ AF 



2 F F 



+- 



UJ 



1 

2F 



, 2 da; 9/ dV „ • / F 
d0 d<p dj> I 2F 



A0 



(2.20) 



$4 + $ 



AF 



(2.21) 



AF AF .„ AF uj ■ ■ 1 /; 2 du; 9/ 
^ + 2tf— + (F-p) — + -0A0 + -U 2 ^ + J 



2^]A0 



(2.22) 
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and 



AF = 6 



$ H + 4H$ H + ~$ H - F$ A - Uh + AH 2 -^)<$> a 



(2.23) 



In the de Sitter background (H ,<p ), these equations assume the considerably simpler 
form 



a 2 2uri 



2uJn 



AR 



AF + 3F AF + \—-4H£\AF + ^AR = 0, 



Ht + SHq Ht H — r i?T = , 

a 2 

1 / AF AF 
-<&h + H $ A — 7: H —— 



1 AF 



H 



$4 + $ 



2 F 
AF 



H 



"0 



(2.24) 
(2.25) 
(2.26) 

(2.27) 
(2.28) 
(2.29) 



$ H + 3H <b H - H^ A - 3H 2 $ A 
to first order, whereas 
AF = 6 



1 AF 



2 F n 



AF 3F 2 AF 

-no H " TT - 



■2 k 2 ( k 2 

$h + 4F $ H + Fo$a + — 



4F n $ A 



(2.30) 



(2.31) 



To first order and in the absence of ordinary matter [35], vector perturbations do not 
have any effect. Expanding de Sitter spaces with H > are always stable, to first order, 
with respect to tensor perturbations [36J, as can be seen from eq. (|2.18j) . On the other 
hand, contracting de Sitter spaces with H < are always unstable {TO- and will not be 
considered further. There remain scalar perturbations, which we set out to examine. 

In modified gravity theories with / = /(F), 0=1, and $rr ^ 0, eqs. ()2.28|) and 
EZB) yield 

AF 

®h = ®A = ~ W , (2-32) 
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whereas eqs. (|2.31|) and (|2.32|) yield 



AR = 6 



a 2 



AHl 



(2.33) 



and a = a e Hot is the scale factor of the unperturbed de Sitter space, with ao a constant. 
In the de Sitter background the gauge-invariant variables reduce, to first order, to 



A0 



AR = 5R , AF = 5F , Af = 5f. 



Since the function F depends only on R one has 

f RR 



where 



and therefore 



F 

fRR 

AR = 



F 

d 2 f 



AR 



dR 2 
' f. 



Ro 



H 



RR 



The perturbations <3># = $,4 then evolve according to eq. ()2.30|) . which becomes 



A- 2 



F 
3/ 



RR 



H 







(2.34) 

(2.35) 

(2.36) 
(2.37) 

(2.38) 



where a = a e Hot , with ao a constant. At late times the term k 2 /a 2 can be neglected 
and stability is achieved if the coefficient of in the last term of the left hand side of 
eq. ()2.38|) is positive or zero: upon use of the value of the Hubble parameter given by 
eq. ()2.5|) for the unperturbed de Sitter space, this condition reduces to 



F 2 



2fol 



RR 



Fol 



> 



RR 



(2.39) 



This inequality was presented in a previous communication without details of the deriva- 
tion [22] ■ We now comment on the physical meaning of the approximation leading to 
f!2.39|) . Information on the spatial dependence of the inhomogeneous scalar perturba- 
tions are contained in the eigenvector k of the spherical harmonics, and the fact that the 
only term containing k in eq. ()2.38|) becomes negligible as time progresses in a de Sitter 
background implies that the spatial dependence effectively disappears from the analysis. 
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One may be tempted to conclude that the stability condition (J2.39)) with respect to in- 
homogeneous perturbations could be obtained in a much quicker way by considering the 
simpler homogeneous perturbations: this would be incorrect as one would not be able to 
guess a priori, in a homogeneous perturbation analysis, the structure of eq. ()2.38|) and 
the fact that the spatial dependence disappears. Furthermore, in the parallel case of 
scalar-tensor gravity, the stability condition with respect to homogeneous perturbations 
differs from the corresponding one for inhomogeneous perturbations, and one would 
expect the same to happen for modified gravity. This is the subject of the next section. 



3 Homogeneous perturbations in modified gravity 
and in scalar— tensor theories 

We now derive the stability condition of de Sitter space with respect to homogeneous 
perturbations in modified gravity, in order to compare it with the result (|2.39J1 of the 
previous section. The field equations reduce to 



H 



H 



1 fRF-f 



3F 



3HF 



±(f-hf) 

2F\ J 



By assuming that H(t) — Hq + 5H(t) and using the first order expansions 



(3.1) 
(3.2) 



R = R + 6R, 5R = 6\6H + 4H SH 
F = F + f RR 5R , f = f + F 5R , (3.3) 

and eq. ()2.5|) . one obtains the evolution equation for the homogeneous perturbation 5H 

fo 



SH 



I / 1 , ) — — I SH + — ( — ^ — } SH 0. 

6H F J 3 \f RR F 



(3.4) 



The ansatz SH = e e yields an algebraic equation for s with roots 



s± = - 
2 



fo 



2H F 



± 



fo 



2H F 



4 (F^_2f 1 
3 \fnR F 



(3.5) 
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Assuming f > and H > 0, if also F > 0, then —fo/ (2H F ) < and there is 
stability if and only if 

2/ , , 

7 ^ - ' ( 3 - 6 ) 

which is equivalent to the condition ([2.39)1 . If ()3.6|) is not satisfied, the root s + is real 
and positive, corresponding to an unstable mode growing exponentially in time. 

The case F < does not correspond to a de Sitter solution when f > because 
eq. ()2.5|) . which reduces to RoF = 2f cannot be satisfied in this case. 

Why the stability condition with respect to homogeneous perturbations coincides 
with the corresponding stability condition with respect to inhomogeneous perturba- 
tions? A naive answer would be that the spatial dependence of the inhomogeneous 
perturbations can safely be eliminated in the analysis of eq. (J2.38)) or, in other words, 
inhomogeneities are redshifted away, and the results must coincide. This would be 
intuitive: even initial anisotropics are known to be smoothed out by de Sitter-like ex- 
pansion [37] in general relativity and in some scalar-tensor cosmologies [3B], but this 
is not the correct explanation: in fact, if it were true, it should hold also in the case 
of scalar-tensor gravity in which / (cf>, R) = ip (</>) R, but this is not the case as we are 
going to show. The stability condition with respect to inhomogeneous perturbations 
in scalar-tensor gravity has been derived in Ref. [SH] by analysing the equation for the 
gauge-independent Bardeen potentials and for the Ellis-Bruni variable A0, 



A0 + 3H A(p + 



n'2 



to V" _|_ 

2 V + 



Fq 



Acp = 



where 



f<f>R 



d 2 f 



d<pdR 



f" 
Jo 



d 2 f 



(<Ao,fio) 



(4>o,Ro) 



This equation is obtained from eqs. (l2~2H) - (l2~3Tl) if 1 + 3/| fl / (2^^) ^ 0. The 
condition that ensues is 



f 

iSL — V" 4- 
2 V + 



Fo 



OJo 1 + 



2^0 Fo J 



< 



For the sake of comparison, let us derive the corresponding stability condition 
spect to homogeneous perturbations in scalar-tensor gravity. Assuming that / 
if)(<j>)R, the homogeneous perturbations 



(3.7) 

(3.8) 
stability 

(3.9) 

with re- 

[</>,R) = 



H(t) = H + 5H(t) , <f>(t) = 0o + H{t) 



(3.10) 
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satisfy the first order evolutions equations 

5h = ~~ Wo ~ Ho ^) ' (3 - U) 

54> + 3H 6j> + (2V " - Q 6<j> = . (3.12) 

ZUJo 

By contrast, in general relativity with a minimally coupled scalar field the perturbations 
have no effect to first order. This is due to the non-canonical form of the effective 
energy-momentum tensor of the scalar field appearing in the left hand side of the field 
equations of scalar-tensor gravity when these are written in the form G a b = ^T^^ [0] 
(see Refs. [33110] for a discussion). The stability condition with respect to homogeneous 
perturbations can be read off of eq. (|3.12j) . 

lo y" 

< . (3.13) 

The stability condition ()3.9|) with respect to inhomogeneous perturbations is more re- 
strictive than 1)3.13)1 and, in spite of having neglected a term k 2 /a 2 = k 2 e~ 2Hot / in 
eq. (|3.7|) . the final stability condition ()3.9)) retains a memory of the spatial dependence of 
the inhomogeneous perturbations, which is instead lost in the homogeneous perturbation 
analysis leading to (|3.13|) . A priori, one should expect a similar situation for modified 
gravity, and the fact that the two stability conditions coincide for these theories appears 
to be coincidental. 

For the particular class of scalar-tensor theories described by the action 



S = d 4 x 



0i?-^^ 6 V a 0V 6 0- V{<j>) 



(3.14) 



and containing a single arbitrary coupling function u ((/)), the stability conditions ()3.9)) 
and ()3.13jl coincide [H]. However, these two conditions fail to coincide in the general 
scalar-tensor theory (jl.5|) because, in the right hand side of eq. (|2.24jl . curvature pertur- 
bations AR act as sources for the perturbations A0 (which, in the de Sitter background, 
coincide with S(p). On the contrary, such a term is absent in the homogeneous pertur- 
bation analysis of the Klein-Gordon equation ()2.4j) . and an analogous term does not 
appear in eq. (|2.38|) for the gauge-independent Bardeen potential in modified gra- 
vity because there is no scalar field in this case and = — eq. (|2.24|) then becomes 
homogeneous. 
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4 Application to specific modified gravity scenarios 



We now proceed to apply the stability condition (|2.39|) to certain specific modified 
gravity scenarios that have been proposed in the literature, for which a de Sitter space 
is relevant. 



R 

This theory d 123 121 CHI EH] > with the mass scale A*o - #o ~ 1(T 33 eV, was the 
first candidate proposed to explain the cosmic acceleration, and it is known to be subject 
to an instability that develops on a time scale of order 10~ 26 s |24j . For this theory, the 
stability condition ()2.39|) reduces to 

1 + (4.1) 

and it is clear that by taking fi w H with i? = 12-ff, 2 , the stability condition is 
impossible to satisfy. More precisely, the condition for the existence of de Sitter solutions 
is R = \/3/i 2 (see also Ref. [2]) and the stability condition ()2.39j) reduces to 8/3 < 0, 
which obviously cannot be satisfied. We stress that this instability of de Sitter space 
arises in the gravitational sector of the theory, while the instability discovered in Ref. [21] 
arises in the matter sector and would disappear in vacuum. The instability of de Sitter 
space can be stabilized by adding a term of the form eR 2 with < e < // 4 to the 
Lagrangian density, as shown in the following. 

• f(R) = R-^ + aR 2 

In this theory [H] the condition for the existence of de Sitter solutions ()2.5j) becomes 

na us! 

Ro = V3 fJl 2 , (4.2) 

independent of the parameter of the quadratic correction — therefore this condition holds 
true also for a = as seen in the previous case. For general values of the parameter a, 
upon use of eq. ()4.2j) . the stability condition (j2.39|) reduces to 

' > (4.3) 



2 



3v3a£* 

and therefore de Sitter space is stable if a > (3a/3/x 2 ) and unstable if a < (3v^3 jj 2 ) 
(in particular for a = 0, which is the previous case). Therefore, adding a quadratic 
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correction with a < (3-\/3/i 2 ) (and in particular with a negative a, which reinforces 
the effect of the term —fi 4 /R) leads to instability. Now, if fi ~ H ~ 10~ 33 eV, the 
parameter a must be larger than ~ 10 65 (eV)~ 2 for stability, which appears to be huge 
in natural units — stability is achieved at the price of fine-tuning the parameters. 



• f(R) = R n 

This modified gravity theory has been pursued in the literature [EJUS!) especially for 
n — — 1 and for n = 3/2, in which case it is is conformally equivalent to Liouville 
field theory [T^]. The model has been used to explain the cosmic acceleration and it 
is also interesting because ordinary inflation with a minimally coupled scalar field and 
an exponential potential (power-law inflation) can be rewritten as a theory f(R) = R n 
[44 . For generic values of n ^ 0, 1/2, 1, the theory yields power-law inflation a oc t a 
with 

-2n 2 + 3n-l 

« = ^ • 4.4 

n — 2 

This theory does not admit de Sitter solutions if n ^ 2 jl?! EH] unless a cosmological 
constant, corresponding to a term with n = 0, is added to the action [47J. However, 
Minkowski space (the trivial de Sitter space) is a solution without cosmological constant 
for any positive value of n [IB]- in fact, the condition ()2.5|) for the existence of de Sitter 
solutions yields nR = 2Rq, which is only satisfied for n = 2 or Rq = 0. The stability 
condition (|2.39|) yields 

2 — n 

Ro > (4.5) 

nyn — 1) 

and is satisfied for 1 < n < 2 and for n < 0; it is identically satisfied for n = 2, without 
imposing constraints on the Hubble parameter H of de Sitter space. The Minkowski 
spaces Rq = are stable for any n > 0. 



• f(R) = R + eR 2 

Quadratic corrections to the Einstein-Hilbert Lagrangian density are motivated by 
renormalizability ^SJ El and higher order corrections are unavoidable near the Planck 
scale jlH] • This theory was used in one of the earliest inflationary scenarios [T7j , not re- 
quiring an inflaton field. The constant has dimensions e ~ M~ 2 , where M ~ 10 12 GeV. 
The condition for the existence of de Sitter solutions allows only the trivial Minkowski 
space H = 0. However, there are non-trivial de Sitter solutions if a cosmological con- 
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stant is added to f(R) [SHI. The stability condition (|2.39|) yields 

1 



e(l + 2eR 



> (4.6) 



which, for Minkowski space, gives stability if e > and instability for e < 0. The case 
e = corresponding to Einstein's theory must be studied separately and it is concluded 
that Minkowski space is stable in this case 

• / (R) = R + eR 2 -2A 

By adding a cosmological constant to the previous theory, de Sitter solutions become 
possible and are given by Rq = 4A, or Hq = a/A/3 as in general relativity, because 
the condition on H does not depend on the parameter e and coincides with the corre- 
sponding condition for e = 0. The stability condition (|2.39j) reduces to e > 0: a positive 
quadratic correction acts in the same direction as the term R in the Lagrangian, whereas 
a negative quadratic correction with |e| arbitrarily small makes de Sitter space unstable. 



• f(R) = R + eR n 

This theory ^3J|5U|, with e ~ M 2<yl ~ n \ where M is a mass scale, comprises quantum 
gravity-motivated corrections to the Einstein-Hilbert Lagrangian for n > 0, as well 
as the theory f(R) = R — fJ> A /R already discussed, or similar theories, if n < 0. The 
condition for the existence of de Sitter space is either Rq = (Minkowski space) or 

K- 1 = (12HJ)- = ^ (4.7) 

for n 2 (the case n = 2 has already been considered). In the case of a non-trivial de 
Sitter space H ^ the stability condition (|2.39|) yields, using eq. (|4.7j) . 

-n 2 + 2n - 2 . t , 

> (4.8) 



en in — 1, 



or en < 0, and therefore de Sitter space is stable if en < and unstable otherwise. In 
particular, it is stable if e > and n < (or if e < and n > 0), which comprises the 
case f(R) = R — fi 4 /R (as n = 2 is excluded from this analysis, these results do not 
contradict the previous statements on the stability of Minkowski space when n = 2). 
Therefore, any theory of the form f(R) = R — / u 2 ^ 1 ~™- ) / R m with m > exhibits the same 
instability in the gravitational sector as the the model f(R) = R — fi 4 /R and, likely, the 
same instability reported in Ref. |23| for the matter sector. 
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f(R)=aln (f) 



This theory [23 E2J EB] , which does not admit a Minkowski space or other solutions with 
vanishing Ricci curvature, admits a de Sitter solution only if R = V2Hq = be b ^ 2 , where 
b is a positive parameter. The stability condition ()2.39|) reduces to 



which, for b > 0, is equivalent to e b < and obviously is never satisfied: the de Sitter 



5 Discussion and conclusions 

The general stability condition (|2.39|) of de Sitter space in modified gravity derived in 
sec. 2 allows one to quickly assess the stability of de Sitter space in specific scenarios. Of 
course, when the number of terms in the Lagrangian density grows, so does the volume 
of parameter space to be searched and this analysis becomes cumbersome — it would 
be greatly helped by an estimate of the range of values of the parameters involved. In 
the future we plan to extend to power-law solutions the study carried out here for de 
Sitter spaces. 

The comparison of the modified gravity results with the analogous results in scalar- 
tensor gravity shows that, although there exists a dynamical equivalence between mod- 
ified gravity and scalar-tensor gravity [S3], this should not be taken too literally. The 
stability conditions with respect to homogeneous and inhomogeneous perturbations co- 
incide in modified gravity but not in scalar-tensor gravity, due to the different detailed 
structure of the equations satisfied by the perturbations. 

The scope of the stability analysis can perhaps be extended to more general theories of 
the form / (i?, R 2 , R a bR ab , R a b c dR abcd ) containing string-motivated corerctions (see, e.g., 
Refs. 1561 IBTj ): there are however doubts on certain choices of the string corrections 
to the Einstein-Hilbert action, due to ghosts or light, long-ranged, gravitational scalars 
that potentially violate Solar System bounds [30] 

Finally, it should be stressed that, while we have analyzed linear stability with respect 
to inhomogeneous perturbations described by gauge-invariant variables, other definitions 
of stability can be considered: for example, stability with respect to black hole nucleation 
[15] or quantum fluctuations [5S] • Sometimes these different definitions yield results that 
are qualitatively similar to our stability condition (|2.39|) (e.g., j3H])- Another possibility 
is to search for a positive-definite energy functional [59], although one should probably 




(4.9) 



space H = ^Jb e 6 / 2 /12 is unstable. 
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look for energies bounded from below rather than positive energies |6Uj ; and so on. 
The fact that these stability criteria are inequivalent is not surprising since the physical 
processes considered are quite different and, even from the mathematical point of view, 
several inequivalent definitions of stability exist for dynamical systems [6*T] . 
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